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We propose a model for the multi-orbital material Sr2CrOs06 (SCOO), an insulator with remark- 
able magnetic properties and the highest T c ~ 725 K among all perovskites with a net moment. We 
derive a new criterion for the Mott transition {U\U2) 1 ^ 2 > 2.5W using slave rotor mean field theory, 
where W is the bandwidth and Uws) are the effective Coulomb interactions on Cr(Os) including 
Hund's coupling. We show that SCOO is a Mott insulator, where the large Cr Uj compensates for 
the small Os U2- The spin sector is described by a frustrated antiferromagnetic Heisenberg model 
that naturally explains the net moment arising from canting and also the observed non-monotonic 
magnetization M(T). We predict characteristic magnetic structure factor peaks that can be probed 
by neutron experiments. 



Transition metal oxides [TJ with partially filled 3d 
and 4d shells have dominated materials research in past 
decades, leading to such spectacular phenomena as high 
T c superconductivity and colossal magnetoresistance. We 
expect the next revolution to involve 5d oxides where the 
interplay of spin-orbit coupling (SOC) and strong corre- 
lations can lead to novel phases of matter. Recent ex- 
perimental efforts have focused on iridium oxides [2H5], 
driven by the possibility of exotic phases (5H5J. Here we 
explore other equally promising experimental systems, 
the 5d oxides containing osmium that have received far 
less attention and provide a rich area for investigation. 

We focus on the double perovskite (DP) family of ox- 
ides with the formula A2BP/O6. A is an alkali or alka- 
line earth metal while B and B' are two transition metals 
(TMs) arranged on a 3D checkerboard lattice. With two 
TM ions, the range of properties span metals to insula- 
tors, ranging from ferromagnets, antifcrromagnets, ferro- 
electrics, multiferroics, to spin liquids |10H13j . 

The most studied DP is Sr 2 FeMoOg, a half metallic 
ferrimagnet with T c « 420K [TTJ [JJ]. It is now well- 
understood that a generalized double exchange mecha- 
nism can explain ferromagnetism with the scale of T c set 
by the kinetic energy of itinerant electrons |15j . From this 
perspective, the observation of an even higher T c &s 725K, 
the highest T c amongst all perovskites with a net mo- 
ment, in an insulator S^CrOsOg (SCOO) [T7] is rather 
puzzling. 

There have been several important density functional 
theory (DFT) calculations of SCOO [T7||25H27], never- 
theless many puzzling questions remain open, (i) Why 
is SCOO an insulator [28]? It is not a band insulator, 
since the bands are partially filled. It is not obviously a 
Slater insulator, given the large moment observed on Cr. 
It is not a priori clear how it can be a Mott insulator 
either, given the weak correlations on Os relative to the 
large bandwidth of 5d orbitals. (ii) Why is there a net 
moment, given that both Cr and Os are in d 3 configu- 
rations? What is the role of SOC on Os? DFT calcula- 
tions 25, 27 find a net moment only if SOC is included, 



but XMCD experiments 17] show insignificant Os orbital 
moment, (iii) What sets the scale for the high magnetic 
T c ~ 725K? DFT does not give a clear answer (25]. (iv) 
Why is the magnetization M(T) a non-monotonic func- 
tion [17j of T? Finite temperature properties are, of 
course, very hard for DFT to address. 

In the Letter we gain insight into all these questions 
within a theoretical framework that describes the hierar- 
chy of energy scales in the charge and spin sectors. 

(1) We use a slave-rotor mean field theory to analyze the 
charge sector of SCOO and show that it is a multi-orbital 
Mott insulator. We derive a new criterion for the Mott 
transition for this system: (U1U2) 1 / 2 > 2.5W where W 
is the bandwidth and Uip) the effective Coulomb repul- 
sion on Cr(Os), which includes the important effects of 
Hund's coupling in addition to Hubbard U. Thus a small 
U2 on Os is compensated by a strong U\ on Cr and drives 
the system into a Mott insulating state. 

(2) We find that the orbital angular momentum on Os is 
quenched, and the Os moment is purely of spin origin, in 
agreement with XMCD measurements [17] . 

(3) We show that the effective spin Hamiltonian for 
SCOO is a frustrated Heisenberg model with antiferro- 
magnetic exchange between nearest neighbor Cr-Os as 
well as next-nearest neighbor Os-Os spins. We find that 
the net magnetic moment M(0) arises, not because of 
SOC, but rather because of a canted magnetic ground 
state, found using both Monte Carlo (MC) simulations 
and a variational approach. 

(4) The scale for the high T c is set by a J\ = 35 meV 
superexchange between Cr-Os near neighbors. 

(5) We find that we can obtain an unusual non-monotonic 
M(T), in agreement with magnetization and neutron ex- 
periments |17j , due to a subtle frustration effect between 
Cr-Os Ji and Os-Os J2. 

(6) Finally, we predict distinct magnetic structure factor 
peaks for our canted ground state. These can be mea- 
sured by neutrons and provide a "smoking gun" test of 
our theory. 

Hamiltonian: The cubic crystal field splits (i-orbitals and 
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FIG. 1: Phase diagram obtained from slave-rotor mean field 
theory, (a) Schematic of Cr and Os orbitals and relevant en- 
ergies, (b) Phase boundary of eq. |2| for the particle-hole 
symmetric case t' = 0, A = Ao, and zi = fio; see text. Com- 
bination of small (72 on Os and large Ui on Cr can drive the 
system Mott insulating. Effect of (c) Os-Os hopping t' and 
(d) Cr-Os charge-transfer energy A. 



only the t 2g orbitals are occupied, with both Cr and Os 
in d 3 configurations. To investigate the metal-insulator 
transition, we consider 

H = -t ^ {4a C jc + M - t' E ( C L C 3« + h.C) 
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-flj^nua + n 2 ia) - A/Xnua ~ n 2 i a ). (1) 

This has orbitally-symmetric onsite Coulomb interac- 
tions that make the problem analytically tractable and 
provide insight, while capturing the essence of strong cor- 
relations in a multi-band system. Hund's coupling will 
be discussed below. Here d\ a (di a ) and c\ a (ci Q ) are the 
creation (annihilation) operators on the i th Cr and Os 
sites respectively. n lia (n 2ia ) are the number operators 
on the i th Cr (Os) site. Note that the label a includes 
both spin and t 2g orbital index, t and t' are the hopping 
amplitude between neighboring Cr-Os and Os-Os sites. 
The effective on-site Coulomb interaction U\(2) is related 
to the atomic charge gap on Cr(Os), as discussed in de- 
tail below. The chemical potential is /x and A, the charge 
transfer gap between Cr and Os. 



Multi-band Mott criterion: We investigate the charge 
sector of (JlJ using the slave-rotor mean field (MF) the- 
ory [30]. This method captures the qualitative physics 
of dynamical mean field theory, and is well suited to 
studying the multi-band Mott transition; it has recently 
been used for various other multi-band system includ- 
ing pyroclores and iron pnictides [3TJ [32] ■ We decom- 
pose the physical electron into an 0(2) rotor that carries 
the charge and a spinon that carries the spin and or- 
bital degrees of freedom. We write d\ t 
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and 



9ia e wnere fj(fi) an d gl(gi) are the fermonic 
spinon creation(annihilation) operators on the i th Cr and 
Os sites respectively. The bosonic rotor operator 8i(<f)i) 
is conjugate to the charge on Cr(Os) at site i. Substitut- 
ing into (JlJ and performing a MF decomposition of the 
hopping terms, we obtain spinon and rotor Hamiltonians 
that we solve self-consistently. 

We calculate (cos 9) and (cos </>) which are used as di- 
agnostics of the metal-insulator transition. When both 
quantities vanish, the system is a Mott insulator with 
suppressed charge fluctuations, while non-zero values 
lead to a metal. For simplicity, we first work at the 
particle-hole (p-h) symmetric point: \i = fio = 3(Ux + 
U 2 )/2, A = A = 3((7i - U 2 )/2 and t' = 0, where we can 
obtain analytical results [33] . We find 



{UxU 2 



a/2 



> 2.W, 



(2) 



where W — 8i is the t 2g bandwidth. This new Mott 
criterion generalizes the well-known result U > W . 

Deviations from the p-h symmetric point, relevant for 
real materials, require numerical solution of the slave ro- 
tor mean field equations. In general, the right-hand side 
of ([2]) is replaced by the the average kinetic energy, which 
depends on the filling and band structure [33]. 

The resulting phase diagram is shown in Fig. [l] We 
see from Fig.JlJc) that increasing t' favors the metallicity 
by introducing an additional route for gaining kinetic en- 
ergy. Fig. [ljd) shows that increasing the deviation from 
the symmetric point A = A — A , for a fixed /i = /io 
and t' = 0, also favors metallicity as it increases charge 
fluctuations within a unit cell [33 . 

Determining the effective U: In order to apply our results 
to SCOO, we must relate the effective interaction Ui( 2 ) 
to material parameters: the Hubbard U , Hund's coupling 
Jh and SOC £. This is done by matching the excitation 
energy in the atomic limit of ([T]) to that of the more 
general atomic Hamiltonian 

H? = {U a ~?,J Ha ) Na{N f l) -2J Ha (S a f 

+ l -J Ha { r Iaf - A a T a -S Q - fi a N a . (3) 

Here the label a = 1, 2 indicates Cr and Os sites respec- 
tively. N is the total number operator, S the total spin 



operator, and T the total orbital angular momentum op- 
erator in the t 2g manifold, The effective SOC A in a given 
TS manifold is a function of T, S and £ [33]. We work 
with a d 3 atomic ground state. 

For the Cr site, we ignore SOC and find the atomic 
excitation gap to obtain Ui = U\ + 2 J hi [33] • The ef- 
fective repulsion is enhanced by Hund's coupling, well 
known (34] [35] for half-filled orbitals. On the Os site, we 
include A2 and find U 2 = U% + 2 Jh2 — 3£ 2 /2 [33], so that 
SOC reduces the atomic charge gap. 

We now can use our Mott criterion for SCOO. For 
Cr, we use U\ s» 5 eV and J hi ~ 1 eV [T]. For Os, 
we use U 2 «,2_eV, J hi ~ 0.35 eV and £2 0.3 eV 
[36] , Thus (IW 2 ) 1/2 ~ 4 eV which is larger than the 
critical value 2.5W » 3.75 eV using the band theory 
estimate W w 1.5 eV [27]. We thus believe that SCOO is 
a Mott insulator. We emphasize that Q is a conservative 
criterion for a Mott insulator for two reasons. First, MF 
theory overestimates the stability of the metallic phase 
(which has long range order in the rotor sector) . Second, 
the slave rotor approach neglects magnetic order, which 
would further help charge localization. 

It is worthwhile to check eq. Q for a ternary perovskite 
NaOs0 3 , where it reduces to U > 2.5W. NaOs0 3 does 
not satisfy this inequality and is thus on the metallic 
side of the transition. This is in agreement with experi- 
ments [22JI37] that show metallic behavior in the param- 
agnetic state, with spin density wave magnetism at low 
T. On the other hand, fo£a DP like SCOO, it is the ge- 
ometric mean ofU± and U2 that determines the^effective 
strength of correlation. Thus the small onsite U% on the 
5c? site Os is compensated by the large U\ on the 3d site 
like Cr to drive the system Mott insulating. 

Effective Spin Hamiltonian: Having shown that charge 
is localized, we turn to the spin sector. For simplicity, 
we take a strong coupling approach, although in reality 
SCOO is likely an intermediate coupling material (see 
Concluding section). Given the large Jh on Cr, the three 
spins align to yield S=3/2. We might expect that a large 
SOC Al • S on Os might split the 5d t 2g orbitals into 
J = 3/2 and J = 1/2 manifolds. However, in the special 
case of a d 3 configuration, the orbital angular momentum 
is quenched. Thus we also get a S=3/2 spin on Os. While 
the absence of an orbital moment on Os is consistent with 
XMCD data [17] . we must now explain the net moment, 
since SOC cannot be responsible for this. 
The effective spin Hamiltonian is given by 

(ij) 

Here S^ 2 ' is a classical S=3/2 spin on the i th Cr(Os) 
site; see Fig[2pt>). Given the general rules of superex- 
change for half-filled orbitals, J x between Cr-Os and J 2 
between Os-Os are both antiferromagnetic. We have re- 
tained next nearest neighbor Os-Os couplings because 
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FIG. 2: (a) Canted ground state obtained from MC simula- 
tions (4096 spins, J2/J1 = 0.4). (b) Schematic showing the 
antiferromagnetic interactions, (c) M(T) for J2/J1 = 0.4. 
By comparing with experimental T c we obtain Ji=35 meV. 
Note that M(T = 0) is finite and M(T) has an unusual non- 
monotonic behavior in agreement with experimental data in 
inset [17] , 



5d orbitals are much more extended than 3d orbitals. 
Ji tends to order the spins into a Neel antiferromagnet, 
while J 2 frustrates this ordering. Note that J 2 by itself 
is also frustrated since Os spins are on an FCC lattice. 
Given the complexity of estimating superexchange micro- 
scopically [38 , we treat them as parameters and study 
the magnetic properties as a function of J 2 / J\. 

Magnetic Structure: In Fig. |2ja), we show the mag- 
netic ground state obtained from MC simulation of 
H c g for J 2 /Ji — 0.4. The magnetic structure fac- 
tor has two large peaks at qo = (7r,7r,7r)/a and qi = 
{it, 7T, 0)/a, where a is the nearest neighbor Cr-Os sep- 
aration. We gain insight into this ground state us- 
ing a variational analysis. Fourier transforming Q we 
get H cS = E 9 -0 f (q) ■ H(q) ■ t/>(q) where H n = 0, 
#12 = #21 = Ji(cos{q x a) + cos(q y a) + cos(g 2 a)), H 22 = 
J 2 (cos((q x - q y )a) + cos((q v - q z )a) + cos({q x - q z )a)) 
and tp^ = (S 1 (q)*, S 2 (q)*). For each q we diagonalize 
-ff(q), and find that qo = (tt, 7t, 7r)/a, the usual antifer- 
romagnetic Neel ordering of Cr and Os spins, minimizes 
the energy. While it satisfies IS^ 2 + |Sf| 2 = 2S 2 , it does 
not satisfy the constraints \S}\ 2 = S 2 and |S 2 | 2 = S 2 for 
■h -I. : 0. 
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TABLE I: Oxides of Cr and Os, where I denotes Insulators, 
M metals, AF antiferromagnets, and Fi ferrimagnets. The 
Weiss constant 6 is not available for materials with high T c . 



Upon generalizing the solution to include two differ- 
ent q's that minimize the energy and also satisfy the 
constraints, we find that qo = (tt, 7r,7r)/a and qi = 
(tt,tt, 0)/a in agreement with the MC results. Our pre- 
diction of the structure factor peaks can be checked by 
neutron scattering, and should prove to be an important 
test of our theory. 

Magnetization M(T): We calculate M(T) in 3D using 
finite temperature classical MC simulations. Fig. [2jc) 
shows that (i) M(T = 0) is finite, (ii) M (T) has an 
unusual non-monotonic function of temperature, (iii) a 
large value of T c for Ji = 35 meV and J2/J1 — 0.4. 
Remarkably, all of these results compare very well with 
experimental data shown in the inset of Fig. [2^c) [T7] . 

The non-zero value of M(0) arises from the canted na- 
ture of the ground state discussed above; see Fig. [2ja). 
There are some quantitative differences with experiment, 
e.g., M(T = 0) R3 0.75/i_B and it peaks at w 0.9fiB while 
the simulation give M(T = 0) « 0.5 fis, which could be 
due to any one of several effects ignored in our minimal 
model, like higher order effects of SOC or distortion from 
the cubic lattice. 

The non-monotonicity of the net M (T) arises from 
a competition between the magnetizations on the two 
sublattices: for essentially antiparallel moments M — 
Mcr — Mos- At low T, the Os moments are stuck in 
the canted state and Mo s does not change significantly 
with temperature. However Mo drops faster leading to a 
decrease in M. This trend continues till about 150K. At 
intermediate temperatures, Os spins have a smaller spin 
stiffness compared to Cr spins because J2 acts only on 
Os sublattice. Thus Os spins become more "floppy" and 



depolarize faster leading to an increase in M with tem- 
perature. Close to T c both Cr and Os moments drop very 
rapidly and the net moment goes to zero at T c . This pro- 
vides a qualitative understanding of the non-monotonic 
behavior. The behavior shown in Fig. |2jc) depends sen- 
sitively on the choice of J2/J1, and M(T) for different 
J2/J1 is discussed in ref. [33] . 

Discussion: For weakly correlated materials, the mag- 
netism is of the Slater/spin-density wave form where the 
state above T c is a metal and the charge gap is determined 
by magnetic ordering, In a local-moment Mott insulator, 
the charge gap scale is set by the Coulomb interaction 
and is much larger than T c , set by superexchange. SCOO 
is in effect just on the Mott side of this crossover between 
Slater and Mott magnetism. By combining a weakly in- 
teracting 5d Os with a strongly interacting 3d Cr, we are 
able to tune the interaction strength to an intermediate 
regime and, therefore, maximize the T c . We already see 
some evidence for the emergence of such a crossover be- 
havior: NaOsOa is a Slater insulator with a T/v sa 410 
K, while LaCrOa a strong coupling Mott insulator with 
a 3.5 eV charge gap and Tjv ~ 320 K. SCCO is at an in- 
termediate coupling on the Mott side, and this generates 
a very high ordering temperature. 

Moreover, the insulating and magnetic mechanisms 
that we have unravelled in SCOO are an important step 
for understanding the trends, such as metal to insulator 
transitions, in materials with fixed Os and varying 3d 
ion (Sc-Cr-Fe-Co), as well as fixed Cr and varying 5d ion 
(W-Re-Os) [see Table 1]. 

Conclusion: We have established a general theoretical 
framework for understanding Mott insulators in DP's 
with half-filled bands. It can be easily generalized to 
other materials like Sr 2 MnOsOg, Sr 2 FeOs06 [39] . An im- 
portant open question is the effect of doping these multi- 
band Mott insulators. Away from half-filling, SOC will 
become important, and the interplay of strong SOC and 
interactions may produce exciting new physics. 
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